Parity- Violating Nuclear Force 
as derived from QCD Sum Rules 



W-Y. P. Hwang and Chih-Yi Wen 
The Leung Research Center for Cosmology and Particle Astrophysics, 

Center for Theoretical Sciences, 
Department of Physics, and Institute of Astrophysics, 
National Taiwan University, 
Taipei, Taiwan 106, ROC 

January 26, 2008 

Abstract 

Parity- violating nuclear force, as may be accessed from parity violation stud- 
ies in nuclear systems, represents an area of nonleptonic weak interactions which 
has been the subject of experimental investigations for several decades. In the 
simple meson-exchange picture, parity-violating nuclear force may be param- 
eterized as arising from exchange of w, p, lj, or other meson(s) with strong 
meson-nucleon coupling at one vertex and weak parity- violating meson-nucleon 
coupling at the other vertex. The QCD sum rule method allows for a fairly 
complicated, but nevertheless straightforward, leading-order loop-contribution 
determination of the various parity-violating MNN couplings starting from 
QCD (with the nontrivial vacuum) and Glashow-Salam- Weinberg electroweak 
theory. We continue our earlier investigation of parity- violating ttNN coupling 
(by Henley, Hwang, and Kisslinger) to other parity-violating couplings. Our 
predictions are in reasonable overall agreement with the results estimated on 
phenomenological grounds, such as in the now classic paper of Desplanques, 
Donoghue, and Holstein (DDH), in the global experimental fit of Adelberger 
and Haxton (AH), or the effective field theory (EFT) thinking of Ramsey- 
Musolf and Page (BP). 



PACS: 24.80.+y (Nuclear tests of fundamental interactions and symmetries), 
24.85. +p (Quarks, gluons, and QCD in nuclei and nuclear processes). 
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1 Introduction 



Parity violation (PV) studies in nuclear systems, such as the asymmetry [TJ [2] in 
p + p — > p + p, the photon circular polarization measurement in n + p — > d + 7, 
and PV studies in 18 F and 2l Ne, offer a means of determining the parity-violating 
nuclear force, which represents a special category of nonleptonic weak interactions 
accessible experimentally. In the simple meson-exchange picture, parity-violating 
nuclear force arises from exchange of it, p, 00, or other meson(s) with strong MNN 
coupling at one vertex and weak parity- violating MNN coupling at the other vertex. 
In 1980's, parity-violating meson-nucleon couplings have been estimated primarily on 
phenomenological grounds, such as in the classic paper of Desplanques, Donoghue, 
and Holstein (DDH) [3J. A global fit to obtain these PV meson-nucleon couplings, 
making use of the various experimental data available at the time, was performed 
also (in 1985) by Adelberger and Haxton (AH) pE]. Nevertheless, progress [5], both 
experimental and theoretical, has been slow since then, although the interest in the 
problem arose to some extent in view of the adopted effective field theory for the few 
nucleon systems. 

In 1980's, the standard model of particle physics, which consists of Glashow- 
Salam- Weinberg [GSW] electroweak theory and quantum chromodynamics [QCD], 
had meanwhile been well established. Although the nonperturbative feature of QCD 
manifests itself in formation of hadron structure and makes it very difficult to pre- 
dict quantitatively hadron properties including both the strong parity-conserving and 
weak parity-violating meson-nucleon couplings, the method of QCD sum rules[6] of- 
fers a systematic way for taking into account effects as caused by the nontrivial nature 
of the QCD ground state or the QCD vacuum. The QCD sum rule method has be- 
come quite popular starting from 1990's (till now). In view of the laborious efforts in 
lattice gauge theory, the QCD sum rule method does offer an alternative avenue to 
obtain the various predictions, to some extent in the spirit of QCD. 

To " complete" the study of nuclear parity violations, we need to investigate both 
the parity-conserving and parity- violating meson-nucleon couplings, using the conven- 
tional scheme of nuclear forces. As already shown in another paper [TJ, it is possible 
to use the method of QCD sum rules in external fields [8] to determine the strong 
tcNN, pNN, and uNN couplings. The purpose of this paper is to present a QCD 
sum rule determination of the parity- violating pNN and uNN couplings, which turns 
out to be fairly complicated but nevertheless straightforward. Our present results, 
together with a previous study of the parity- violating ttNN coupling [TU], allow for 
a direct connection between parity-violating nuclear force and the standard model, 
despite the complication of QCD. 



The external- field QCD sum rule method [8 J have been used to treat the strong 
Ti NN coupling [TT] . the weak parity- violating nNN coupling [TP] , and the strong 
pNN and uNN couplings [7J. In all cases considered, quantitative successes have been 
achieved mainly because the nonperturbative effects of QCD, as expressed in terms 
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of induced condensates, have been taken into account and are found to be of critical 
importance. Our efforts to treat hyperon weak decays remain in progress |12j . with 
fairly encouraging results (which are beyond the scope of the present article). There- 
fore, it is natural to follow the same method to treat the weak parity-violating pNN 
and ujNN couplings. Nevertheless, it remains difficult to assess properly whether 
our present calculation may also succeed in a quantitative manner, since after all the 
nonleptonic amplitudes in question are "notoriously difficult to calculate" (quoting 
the phrase from the referee on the early version of one of my early articles). However, 
we have good reasons, based on successful experiences on the method of QCD sum 
rules in general, to believe that the complicated task carried out in this research could 
be the first important step in establishing a benchmark in the (future) quantitative 
treatments of this difficult problem. Just like many other nonleptonic weak decays 
where relevant data exist, we already have experimental information [U [5] which we 
may use to test our theoretical predictions. 

For the sake of completeness, we begin by outlining a few ingredients regarding 
the external-field QCD sum rule method, without detailed qualifying statements. The 
external vector field is expressed as 

Zfj, = — -Z^ u x v . (1) 

We attempt to determine the following polarization function for a nucleon in a small 
(classical) external vector- meson field Z^ (where Z represents either a p or uj meson): 

i J d A x e^ x (0\T VN (x)f} N (0)\0) z = U(q) + hZ^W^q), (2) 

where the numerical constant h is the coupling between the external field and the 
up (u) quark field; more explicitly, the coupling between the external field and the 
down (d) quark is —h for p meson and +h for uj meson. The standard form for the 
composite operator r]^{x) is adopted [9J: 

Vp (x) = e abc [u aT (x)C lfl u b (x)} l5 Yd c (x), (3) 



Vn {x) = e abc [d aT (x)C llM d b (x)} l5 ru c (x), (4) 

which transform like the proton and neutron fields, respectively. Here u a (x) and d a (x) 
are the up and down quark fields with the superscript a the color index and C is the 
charge conjugation operator. 

At the hadronic level, we define the parity- violating meson-nucleon couplings re- 
lated to p and uj mesons in the standard manner [3]. 



+^iV(r x tff—JL^N - N[hlr, + h'y^YlsN. (5) 
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Here the superscripts in the couplings h° p , h p , h 2 , h' p l , h^, and refer to the isospin 
character as the weak interactions do not observe isospin symmetry. We may use the 
above expression to write down the polarization function 11(g), with the paricular 
piece which is proportional to h, at the hadronic level. It is clear that we need to 
combine different channels (i.e. p — > p, n — > n, p — > n, and n — > p) in order to extract 
these PV coupling constants. 

At the quark level, we use weak interactions as described by the Glashow-Weinberg- 
Salam (GSW) electroweak theory to determine the polarization function, obtaining 
in general three-loop diagrams which require regularizations. Here we adopt dimen- 
sional regularization in the minimum-subtraction (MS) scheme and introduce suitable 
counter terms in defining the renormalized operators [r] p ]R and [^„]r. As mentioned 
above, we may parametrize the same polarization function phenomenologically, mak- 
ing use of the parity-violating pNN and uNN interaction Lagrangians. Comparing 
the results obtained through the two ways of evaluating the polarization functions 
(i.e., at the quark level using GSW theory and QCD and at the hadronic level involv- 
ing pNN and uNN couplings), we have a definitive way in extracting the weak PV 
pNN and uNN couplings — the primary objective of this paper. 



2 The parity-violating pNN and ujNN couplings 

As a special feature in relation to the method of QCD sum rules, the kinematic 
variable q 2 of the polarization function n(g), as translated into the choice of the Borel 
mass squared M 2 , which is in the vicinity of slightly above 1 GeV 2 . Such choice of the 
Borel mass is to ensure the approximate validity of the operator-product expansion 
(OPE) augmented with power corrections (as due to the various condensates). In 
other words, perturbative QCD corrections to the coefficients in such OPE are in 
principle there but are presumably suppressed by choice of the Borel mass M. Unlike 
what has been involved in most phenomenological approaches to the problem where 
some effective weak Hamiltonian at the energy scale relevant to the hadron must be 
directly invoked, we have in the QCD sum rule method the nice feature that the GSW 
electroweak theory is called for at the scale set by the Borel mass squared M 2 where 
effects to order 0(Gpas) are suppressed (due to the running of the strong coupling 
as) and it is the intrinsic smooth extrapolation of the results to lower q 2 which helps 
to explain the successes of the predictions. 

Accordingly, it should be possible to use, in the QCD sum rule method, the 
effective GSW electroweak lagrangian at tree level while leaving terms in 0(Gpcts) 
as corrections. That is, we may use, as a good starting point, 

Gp 

L W eak = --^[wy^l - 7 5 )ddy(l - 7 5 )w + 

9i7^(An - B qi -f 5 ) qi q 2 Y{A q2 - B q2 ^ 5 )q 2 ], (6) 

qi,q2=u,d 
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where Gf is the Fermi coupling constant (Gf = 1-166 x 10~ 5 GeV~ 2 ) and the other 
constants are defined through A u = | — |sin 2 ^ , B u = | , A d = — | + | sin 2 ^ 
and I? ' = — | with sin 2 9w the electroweak mixing parameter (sin 2 9w = 0.2315). 
Note that the first term in the above equation comes from exchange of the W ± boson 
while the second term is due to Z Q exchange. In view of our primary goal which is to 
identify the sizable role played by the various condensates, it is clear that the issue 
regarding the renormalization of the tree lagrangian (as would be suppressed by the 
choice of the Borel mass) is of secondary importance and could be addressed when 
higher precision (in theoretical prediction) is called for. As shall be explained later, it 
is in fact sufficient to even use the current-current form since the difference from what 
we may obtain by employing the renormalizable gauge throughout, such as the 
gauge, does not matter at the end as far as our final QCD sum rules are concerned. 



2.1 Renormalization of composite operators 

Inclusion of the weak interaction in the polarization function leads to sub-divergences 
even in the lowest order 0(Gp)- Such divergences may be removed (or regularized) by 
making use of a suitably defined renormalized composite operators [77] r and [77] Such 
renormalized operators may be obtained by considering the four-fermion interaction: 



J d d x d d y e^ x+ ^ (</(0)/(0K(x)4(y)) = ^^p.,^,. + 

%8**(d-4) 7? [ip + q)p{P + ^ + ^ + VI X 

{\YlM U ~ B u l5 )p] t , t [rY(A d - B d l5 )q] fj + 

[7 P 7 M (1 - 75)p] i 'j7 <T 7' 1 (l - 75)5]/^} + finite parts, (7) 

where the first term is the free piece while the second term, which is divergent at 
d — 4, comes from the insertion of the weak interaction. This equation enables us to 
define the renormalized operator [u(x)d(x)]n. Using the minimum-subtraction (MS) 
scheme, we obtain the renormalized operator [u(x)d(x)]n as follows: 



<{x)d){x) = WH^hr) ■ V 



R ' v 'J 

{[YlM U ~ B u l5 )u a (x)Urr(A d - B d l5 )d b (x)] 3 + 
[7^(1 - 75 K(:r)M7'V t (l - 75)^)],}. (8) 

Here and what as follows, the notations are defined in accord with the previous 
equation such as (d p d a + \g P a^) as coming from \(p + q) p (p + q) a + \{p + q) 2 g pa \. 
Similarly, the renormalized operators [u(x)u(x)]r and [d(x)d(x)]R are given by 
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Ui(x)u b Jx) 




(9) 



d*(x)d b Ax) 




(10) 



The renormalized composite operators [i] n ]ji and [i] p ]r may now be deduced di- 
rectly from the above three equations. This corresponds to introduction of the counter 
terms which cancel the sub divergences. We find that [t)]r is given by 



An analogous expression for [rj n ]R may be obtained by exchanging all the w's and gTs. 
We note that, once the sub divergences are removed in this way, the overall divergence 
(associated with the polarization function) becomes local and thus is removed upon 
Borel transform (as employed in the context of QCD sum rules). The difference our 
result from what we may obtain by working throughout with the renormalizable 
gauge also disappears upon Borel transformation. 

2.2 QCD sum rules for the parity-violating couplings 

The various diagrams which we need to consider are illustrated in Figs. 1, where 
the propagator with a thick external-line mark is to be understood as the quark 
propagator in the presence of an external vector field. [3 [8] At the quark level, we 
calculate the polarization function by considering only the leading two terms in the 
operator product expansion — the first one is associated with the operator < 1 > 
and the second with the operator < qa^q >. As explained immediately below, 
the resultant expressions are already extremely complicated algebraically, preventing 



[Vp\r 



e^[u" T (x)Cj»u b (x)) l5 ^d%x) + -J 24 ; 2(rf _ 4) x 

+ i^(D) tf ) x {-[u aT (y)Cn,Yu b (y)} l5l ,d c (x) 

+ {u aT {x)C 1 ,Ylu{A u - B u l5 )u\y)] l5l ,YY(A d - B d lb )d c {y) 
+ [u aT (x)C llM Ylu(l - j 5 )u b (y)] l5llM ri u (l - l 5 )d c (y)}\ y=x . 




(11) 
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Figure 1: Loop diagrams for PV-coupling sum- rule calculations. 



us from carrying out a better and more complete calculation (unless a considerable 
amount of time is invested). 

To get a better feeling towards the extensiveness of the problem, we may look 
at the following expression which is associated with the coefficient of the operator 
< 1 >: 

r d d p! d d p 2 d d p 3 plp%(g - p 3 ) c (p 3 - Pi) d (P3 - P2Y 

J {2ixY{27rY{2iiYplpl{q-p,Y{p l -p,Y{p 2 -p,Y x 

(12) 

which involves taking the trace of the product of 13 gamma matrices, a task that 
needs a good algebraic software package (such as Mathematica). It also involves 
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three integration variables pi,p2 and p% and dimensional regularization is required 
during the integration. To work out the problem, we have chosen to simplify these 
three-loop expressions by devising algebraic programs making use of Mathematica. 
It still takes up considerable amount of computer time in solving the problem. In 
practice, we perform the calculation in two steps: First, we do the integrations in the 
order of pi and p 2 , then p 3 . This needs a liitle program written in Mathematica to 
handle it (using the various formulae suitable in d— dimensions) . The result contains 
about a hundred terms with different tensor structures. Second, we contract the result 
with the gamma matrices and simplifying it. For this step we use the Mathematica 
package "FeynCalc 1.0" (written by Rolf Mertig). The final result is 



-iHv + - 1*Q P + l P Qa)[(A d A u + B d B u ) 

—(A u B d + x + ^L-^ 

[-275 + 108X1 + — - — [4345 - 27tt 2 - 2475X + 486X 2 ]}, 
1 J 165888 1 u 



(13) 



where X = 7 — log(47r) + log(— q 2 ) . The presence of the pole term Xj (d — 4) indicates 
that this diagram contains sub divergences. The sub divergences come from the p\ 
and p 2 integrals. The removal of these sub divergences is done by making use of the 
renormalized operators [i] p ]r and [i] n ] R . 

We have examined the question of how to define 75 in d— dimensions, since the 
definition of 75 is a tricky issue in the d— dimension ('t Hooft & Veltman in 1972). In 
the present case, however, what we may do is to ignore this fact and anti-commute 
it with all the 7^'s to reach the utmost right position. This is what we have done 
in the above expresion, i.e., by moving all the 75 matrix to the right of all gamma 
matrices. After removing all the sub divergences, the difference between using the 
proper procedure and using the naive method turns out to be only a polynomial of 
q 2 . The difference does not contain nonlocal terms like ln(— q 2 ) or 1/q 2 which we must 
handle with care. 

At the hadron level, there are four coupling constants for the p-meson and two 
for the cu-meson so that the calculation should be carried out for the various channels 
including p + np, p~pn, p°pp, p°nn, upp, and unn and suitable linear combinations 
allow for the determination of the various PV couplings. 

We choose to focus on the antisymmetric part of the sum rules (proportional to 
Zfi v with = —Zfi v ), which already contain enough information to determine the 
meson-nucleon couplings. Different tensor structures leads to different sum rules. 
After suitably adding and subtracting between these pNN (uNN) sum rules we 
obtain the following results (in the MS scheme). 
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Ajjg G F /i g 4 



(g 2 - m 2 ) 2 256vr 6 



x 



(31 ln(-g) - 6 ln 2 (-^))(i + sin 2 M (14) 



\%hl G F h g 4 

X 



(g 2 - m 2 ) 2 "~ 256tt 6 

^(59 ln(-^) - 6 ln 2 (-^))(i sin 2 M (15) 

fcj = (16) 
W G^/i ._ . g 2 



X (99) 7^-7 x 



2m(g 2 - m 2 ) 2 ^2 432vr 4 



(131n(-^)-31n 2 (-^))(^sin 2 M (17) 



/i 2 /i 2 3 

\ 2 N hl G F h g 4 



x 



(g 2 - m 2 ) 2 256tt 6 

(31 ln(-^) - 6 ln 2 (-^))(i + sin 2 W ) (18) 



\%hl G F h g 4 

X 



(g 2 - m 2 ) 2 v/2 256tt 6 

^(35 ln(- - 6 ln 2 (-^))(^ sin 2 M (19) 



where ln(l//Z 2 ) = 7 + ln(l/47r// 2 ). The entity \x is the renormalization scale used in 
dimensional regularization. sin 2 9w is the electroweak mixing parameter in the GSW 
electroweak theory. 

Performing Borel transformation on both sides, taking into account anomalous 
dimensions for the various terms in the OPE, and making use of the continuum 
approximation for contributions from higher exited states, we find 



\ 2 N K ^ Gph M 6 4 

f-e 'm 7 = L 9 x 

M 2 v/2 256tt 6 

W W M 1 

[-62£ 2 (-) + 24F 2 (-, -)](- + sin 2 9 W ) (20) 



\ 2 N h\ _ m 2 G F h M 6 r 4 

1 e = — L~9 x 

M 2 256tt 6 



1 W W M 1 

- [ -l 1 8 B2 (-) + 24F 2 (-,-)](-sin^ w) (21) 
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K = (22) 

O A 2 N _rr? G F h M 4 16 

e ^ = — j^x (qq) t^zi l 27 x 



2m M 2 y/2 432tt 4 

W 

X 2 N h° w _n£ G F h M 



W W M 1 

13^i(^)+6F 1 (^,-)](-sin 2 M (23) 



e = = -L s x 

M 2 v/2 256tt 6 



iy W M 1 

[-62£ 2 (^) + 24F 2 (^, -)](- + sin 2 6 W ) (24) 



Ai/ii G F h M 6 



e a7^ = = -L s x 

M 2 v/2 256vr 6 

5 W M 1 o 

— [-70E 2 (— ) + 24F 2 ( — , — - sin 2 ^) (25) 
36 L K M J A M' ft ;JV 3 ; V ; 

The function L (= ln(M/A QC£ ,)/ln(/i/AQ CjD ) with A QCD = 100 MeV and /2 = 
0.5 GeV) is introduced to take care of the anomalous dimensions. Unlike our analysis 
of the QCD sum rules for the strong couplings g p and g w [7\, loop integrations in the 
weak parity- violating case yields sub divergences in the form of p 2n [ln(— p 2 //2 2 )] 2 with n 
some integer which, upon Borel transform, yield terms proportional to M 2n ~ 2 ln(M 2 /p, 2 ) 
in the QCD sum rules. Thus, the value for the dimensional parameter Jx is of some 
numerical importance. The standard choice p, = 0.5 GeV is used in this paper. 

Furthermore, the quantity W (with the standard choice of W = lAbGeV) is the 
threshold used in the continuum approximation, in which the contributions due to 
the excited states and the continuum are approximated by what may be obtained at 
the quark level (by use of QCD). The continuum approximation introduces, into the 
sum rules (20)-(25), the functions Ei and Fj which are defined in the Appendix. 

We should also mention that the QCD sum rules which we have obtained for PV 
couplings are basically those of leading order in QCD (with the nontrivial vacuum 
structure). It is a tremendous task to try to include a sufficient number of terms which 
involves the condensates of higher dimension in nature. Nevertheless, it is important 
to emphasize that the QCD sum rule approach is a definitive (deductive) procedure for 
evaluating the various diagrams based on QCD, contrary to the qualitative or semi- 
quantitative nature of the earlier DDH approach [3]. In other words, the present 
QCD sum approach can be improved upon, albeit fairly complicated, order in order 
in QCD and with increasing dimensions. 



3 Numerical analysis 

The input parameters for our numerical analysis are those commonly adopted in 
standard QCD sum rule analyses [6J l8l [TO]. 
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X 2 N = 1.2 x 1(T 3 Gel/ 6 , a = 0.546 Gel/ 3 , b = 0.47 Gel/ 4 

X = -6 GeV~ 2 , m 2 = 0.8 Gel/ 2 

h = 4.65, # p = 2.79, g u = 8.37 

Variations in some of these parameter values may result in errors of numerical impor- 
tance. Neglect of higher dimensional terms brings in uncertainties which are difficult 
to quantify (unless some such terms can explicitly be taken into account). Neverthe- 
less, we may use our experience from analyzing other QCD sum rules, such as those 
for the strong pNN and uNN couplings |7J, to assess these errors. 

Before performing the standard numerical analysis, we wish to first use the QCD 
sum rules (20)- (25), together with some data from the nuclear parity violation ex- 
periments, to provide some estimates on the various parity-violating (PV) couplings. 
We note that both the couplings h° p and receive contributions from both W ± and 
Z° exchanges while other PV couplings are dictated by Z° exchange. As a result, the 
couplings h° p and are much larger than the other PV coupling constants, which are 
suppressed by a factor (1/3) sin 2 6 w = 0.07. This means that, as an approximation, 
we may neglect the couplings h p , h 2 , h' p and h^. Another important approximate 
relation is h° n = h ,. 

There are many parity violation measurements in processes involving complex 
nuclei, leading to the determination of the quantity X^, which characterizes the 
strength of PV interaction. This quantity, with an experimental value of about 3 x 
10~ 6 [13j, can be expressed in terms of the PV coupling constants: 



X£ = 5.5f w - 0.25g p hl - 0.62g p h° p - 0.05^ - 0.17&X - 0.19#X (26) 

The PV itNN coupling f n may be obtained via the QCD sum rule method and it is 
about (3.0 ±0.5) x 10~ 7 [HJ CED]- Thus, we have h° p « h° w -4.1 x 10~ 7 , by assuming 
that the other PV couplings (h p , h 2 p , h' p , h^) are considerably smaller in comparison. 

Analogously, we may also use the expression of the asymmetry A pp , as observed 
in polarized proton-proton scattering, to determine the PV couplings: 

h 2 

A pp (15MeV) = OMgM + hi) + 0Mg p (h° p + h l p + -|), (27) 

with an experimental value of —(1.7 ± 0.85) x 10~ 7 [3J. It yields a value of about 
— (10 ± 5) x 10~ 7 for h? and h^, the magnitude of which is somewhat larger than the 
previous value but remains consistent in light of large errors. 

We proceed to analyze the various QCD sum rules as a function of the Borel mass 
squared M 2 . The choice of the range for the Borel mass is guided by the study of 
the QCD sum rules for other properties of the nucleon, such as the nucleon mass, 
the magnetic moments, the axial couplings, or the strong and weak ttNN couplings. 
Specifically, it has been found in general that the nucleon sum rules should give rise 
to reasonable predictions in the Borel mass range of 0.9 GeV 2 < M 2 < 1.1 GeV 2 . 
Such choice of the Borel mass squared also enable us to estimate the errors of our 
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QCD sum rule predictions by inferring from the analysis of the other QCD sum rules 
for the nucleon. The error could be as large as about 25 % in some special cases. 



0.080 




Figure 2: The PV coupling h° (or hP u ) is plotted as a function of the Borel mass 
squared M 2 in the range of (0.9 — 1.2) GeV 2 . 



?:K . — i . . r 




Figure 3: The PV coupling h} p is plotted as a function of the Borel mass squared M 2 
in the range of (0.9 - 1.2) GeV 2 . 

In Figs. 2-5, the various couplings are shown as a function of the Borel mass 
squared M 2 (in units of GeV 2 ). It is of some importance to note that the scales in 
these figures are in fact different, resulting in errors of different magnitudes. As long 
as we have faith in these sum rules when a sufficient number of higher dimensional 
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terms are included, it makes sense to make (semi-quantitative) predictions with the 
aid of only a couple of leading terms. The following predictions have been obtained 
in this way: 



K 
K 
K 

v 

p 

K 



h 



-(6.9 ±3.6) x 10~ 7 
-(0.087 ±0.004) x 10" 


(0.23 ±0.02) x 10~ 7 
-(6.9 ±3.6) x 10~ 7 
-(0.03 ±0.04) x 10~ 7 . 



The large errors, such as those associated with hP and , are caused primarily by the 
rapid variation of the prediction in the quoted Borel mass range. The experience with 
the QCD sum rule method, such as in the weak parity- violating coupling f w NN [TO] , 
suggests that inclusion of the next couple of terms of higher dimensions could help to 
smoothen the rapid-varying behavior of the leading term while yielding predictions 
in the same ballpark. 
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Figure 4: The PV coupling h'^ is plotted as a function of the Borel mass squared M 2 
in the range of (0.9 - 1.2) GeV 2 . 

We note that the predicted values for h° p and h^, are compatible with the above 
results as extracted from the experimental values of and A pp {lh MeV) (together 
with the observation from the QCD sum rules that only h° p and are dominant). 

In Table I, we compare the various predictions on parity-violating meson-nucleon 
couplings, including the estimates given by Desplanques, Donoghue, and Holstein 
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Figure 5: The PV coupling is plotted as a function of the Borel mass squared M 2 
in the range of (0.9 - 1.2) GeV 2 . 





DDH (best estimate) 


AH (best fit) 


QCD sum rules 


u 


-> 11(4.6) 


0^ 11(2.1) 


3.0 ±0.5* 




-31 -> 11(— 11.4) 


-31 -> 11(— 5.8) 


-(6.9 ±3.7) 




-0.4 -> 0(-0.19) 


-0.5 -> 0.4(-0.22) 


-(0.087 ±0.004) 


ti 2 


-7.6 -> -ll(-9.5) 


-6.3 -> -10(-7.1) 





h n 
a p 








0.23 ±0.02 


hi 


-10 -> 5.7(-1.9) 


-12 -> 2.6(-5.0) 


-(6.9 ±3.7) 


hi 


-0.8 -> — 1.9(— 1.1) 


-3.1 -> -l.l(-2.4) 


-(0.03 ±0.04) 



Table 1: PV couplings in units of 10 7 . The numbers in the brackets are the "best" 
values. * The prediction is obtained from [14J. 

(DDH) [3], the values by an overall fit to the existing experiments (AH) [1], and 
ours making use of the QCD sum rules. Our results, if taken seriously, sharpen the 
allowed ranges for most of the six couplings. The overall agreement is good, although 
we should note that, from the leading diagrams which we have included, the isotensor 
coupling h 2 vanishes (while it fairly sizable in [31 H] ) and h'^ is small but is different 
from zero (as in [HH H]). [Both h 2 and h'p do not contribute in any major way to the 
existing nuclear parity violation observables. See, for example, the expressions for 
and A pp (15MeV).] 

Owing to the tremendous complications to include a sufficient number of higher 
dimensional terms involving the various condensates, we have relied on the experience 
of analyzing the other QCD sum rules for the nucleon in order to assess the uncer- 
tainties. While it is clearly desirable to improve on these derivations (by including 
more higher order terms, especially for h° p and it is indeed gratifying to note that 
our overall predictions are fairly consistent with earlier results [21 H]. 
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4 Summary 



The problem of parity-violating nuclear force has been a very difficult one. Progresses 
have been made in the past [31 H] in the simple meson-exchange picture, in which 
parity- violating nuclear force arises from exchange of tt, p, u, or other meson with 
strong MNN coupling at one vertex and weak parity-violating MNN coupling at 
the other vertex. The QCD sum rule method allows for a determination of the 
various parity-violating MNN couplings starting from QCD (with the nontrivial 
vacuum) and Glashow-Salam- Weinberg (GSW) electroweak theory. I believe that 
the QCD sum rule methods, which make explicit connection between the underlying 
theory (i.e. the QCD and GSW electroweak theory) and the predictions, offer us a 
systematic method to tackle. We wish to note that our QCD sum rule predictions 
are fairly consistent with earlier results [3j HJ [5]. 
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Appendix: Borel Transformation 

For the invariant functions f(p 2 ) appearing in the polarization function, we may apply 
the following dispersion relation, 

71 JO p 1 + S 

with a necessary number of subtractions. 

Borel transformation is a general tool employed in the method of QCD sum rules. 
It is defined by 

s™ +1 d 

Bf(s)= lim _(-_)»/(*) (29) 
n, s — > oo n - a s 

s/n = M 2 

Apply Borel transformation to the dispersion relation we obtain 

B f(s) = - / e- p /M Im f{p 2 ) dp 2 . (30) 

71 JO 

In general, we may take into account the contribution due to exited states and 
assume that it is equal to the quark-level contribution of f(p 2 ) starting from some 
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cutoff value p 2 = W 2 . That is, we use a modified Borel transformation B w with the 
definition: 

B w f(s) = - / e-v 2 ' M2 Im f(p 2 ) dp 2 , (31) 

7T JO 

which involves some cutoff W. 

We list all the modified Borel transformation functions which we need in the text. 



„ , s 
= B w In — 

jJL 

= _M 2 (l-e- w2 ' M2 ) (32) 



-M 2 E (-) s B^ln 



M«[l- e -^(l + ^)] 03) 



-2M 6 E 2 {^-) = B w s 2 \n^ 



-2M 6 [1 - e ~ w2 ' M2 {l + ^ + -^)] (34) 
1 M 2 2M 4;J V ; 



2M>F l{ ^,-) - S»-m 2 4 
M fi jj, 2 



= 2M*{1 - 7 - e-^l"' + 

+ [l-e-W/-(l + ^)]h^} (35) 

M fi fj, 2 

= _4M 6 {^- 7 -^ 2 / M2 (l + -^-)+^(-^) 
X 2 [ 2 V 3M 2J V M 2; 

Here we have adopted the conventional symbols E ,Ei and E 2 , and introduced two 
new symbols Fi and F 2 . The exponential integral function Ei(x) is defined by 



/CO p t 
—dt (37) 



Ei(x / 

-X 
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